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GRACE/SUSY-loop is a program package for the automatic calculation of the MSSM 
amplitudes in one-loop order. We present features of GRACE/SUSY-loop, processes cal- 
culated using GRACE/SUSY-loop and an extension of the non-linear gauge formalism 
applied to GRACE/SUSY-loop. 

1 Introduction 

Despite its compactness and success in describing known experimental data available up 
to now, the standard model (SM) is considered to be an effective theory valid only at the 
presently accessible energies on account of theoretical problems. Supersymmetric (SUSY) 
theory, which predicts the existence of a partner to every particle of the SM that differs in 
spin by one half, is believed to be an attractive candidate for the theory beyond the SM 
(BSM). The minimal supersymmetric extension of the SM (MSSM) remains consistent with 
all known high-precision experiments at a level comparable to the SM. One of the most 
important aim of the particle experiments at sub-TeV-region and TeV-region energies is to 
probe evidence of the BSM, so search for SUSY particles plays crucial role in it. 

Experiments at present and future accelerators, the Large Hadron Collider (LHC) and 
the International Linear Collider (ILC), are expected to discover SUSY particles and provide 
accurate data on them. In particular, experiments at the ILC offer high-precision determina- 
tion of SUSY parameters via e~e + -annihilation processes. Since the theoretical predictions 
with the similarly high accuracy is required for us to extract important physical results from 
the data, we have to include at least one-loop contributions in perturbative calculations of 
amplitudes. 

Among SUSY particles, only the lightest one (LSP) is stable if i?-parity between usual 
particles and their SUSY partners is conserved, and the others decay without exception. 
Then decay processes should be analyzed precisely in experiments at the LHC and the ILC. 
Recently, we have calculated the radiative corrections to production processes and decay 
processes of SUSY particles in the framework of the MSSM using GRACE/SUSY-loop [TJ [2J 
El IH [5]. In this paper, we show features of GRACE/SUSY-loop, processes calculated using 
GRACE/SUSY-loop and an extension of the non-linear gauge (NLG) formalism [rjl 171 151 M fTTJl 
[TT] applied to GRACE/SUSY-loop [12]. 

2 Features of GRACE/SUSY-loop 

For many-body final states, each production process or decay process is described by a 
large number of Feynman diagrams even in tree-level order. There are still more Feynman 
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diagrams in one-loop order even for two-body final states. For this reason, we have developed 
the GRACE system [TJ] , which enables us to calculate amplitudes automatically. Figure [T] 
shows the system flow of GRACE generically. A program package called GRACE/SUSY-loop 
is the version of the GRACE system for the calculation of the MSSM amplitudes in one-loop 
order, which includes the model files of the MSSM and the loop library. There exist other 
program packages developed by other groups independently for the calculation of the MSSM 
amplitudes in one-loop order, SloopS [H] and FeynArt/Calc [T5] . 



2.1 Renormalization scheme 

As explained in Q], the renormalization scheme adopted for the electroweak (EW) sector 
in GRACE/SUSY-loop is a variation of the on-mass-shell scheme PH [T7J Q21 HB HO] , which 
is a MSSM extension of the scheme in the SM used in GRACE-loop [13]. We impose the 
on- mass-shell condition on gauge bosons (W, Z, 7), all fermions (/), sfermions (/), CP odd 
Higgs {A ), the heavier CP even Higgs {H°), both charginos (xt > xt)i an d the lightest 
neutralino (x?), then these particles have no mass correction in one-loop order of the EW 
sector. There are some freedom in the renormalization scheme of the sfermion sector. They 
are distinguished by different choice of residue conditions, decoupling conditions on the tran- 
sition terms between the lighter and the heavier sfermions, and left-handed SU(2) relations 
in one-loop order. Recently, our calculations have been performed with the scheme in which 
the residue conditions are imposed on all sfermions ^xcept for the heavier stop and sbottom 
(£2 ; ^2)- Corrections of the external line for and &2 become non-zero in this scheme. 

The renormalization schemes adotptcd for the QCD sector are separate ways between 
light and massive particles. Light quarks in the first and second generation and gluon are 
treated in the DR scheme [2TJ[22] as m the convcnsional pcrturbativc QCD. Massive quarks 
in the third generation and gluino are handled by the on-mass-shcll scheme as in the EW 
sector. For the regularization of infrared divergences, the fictitious mass of gluon A is used 
in the previous version of GRACE/SUSY-loop [T]. We have developed a new version of the 
system in which mass-singularities are regularized by the dimensional method [3]. In order 
to refer the ultraviolet and the infrared divergences, we define the notations Cuv = l/ e an d 
Cir = 1/ e, where the dimension of the space-time d is related to e and e as d = 4— 2e = 4+2 e. 



2.2 Non-linear gauge formalism 

In GRACE/SUSY-loop, we use the technique of the NLG formalism in order to confirm the 
validity of calculations by imposing the NLG invariance on physical results. The NLG 
formalism is an extension of the linear i?^-gauge. The gauge fixing lagrangian for the EW 
interactions in the linear i?^-gauge is as follows: 

£gf = -T^W-^^z) 2 -^^) 2 , (1) 



7 



F w ± = d^W ±f "±i^ w M w G ± , (2) 
F z = d^+tzMzG , (3) 
F 7 = d^, (4) 

where G ± and G° stand for the Goldstone bosons which correspond to gauge bosons 
and Z ', respectively. 
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Figure 1: The system flow of GRACE 
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The following NLG functions are introduced to the EW sector of the MSSM Lagrangian, 

F w ± = (d^±ie&A t ,±igcwPZ^W ± ' l ±^ w ^(v + 5 H H + S h h ±ikG°)G ± , (5) 

F z = d^ + ^ z 9 -^-{v + ~e H H + l h h )G , (6) 
F 1 = d»A», (7) 

where v = \J v\ + v\ , M w = gv/2 , M z = gzv/2 , and h° stands for the lighter CP 
even Higgs. They contain seven independent NLG-parametcrs, (a, j3, Sh , Sh, H, £h, We 
perform the numerical tests by varying these parameters. 



2.3 Tests of numerical results 

Since the GRACE system provides numerical results automatically, we have to test validity of 
the results. We adopt four tests for the EW sector and three tests for the QCD sector as in 
Table EJ 





Tests 


Variables 


EW: 


NLG invariancc 






Cancellation of ultraviolet divergence 






Cancellation of infrared divergence 


A (fictitious mass of photon) 




Independence of the cutoff energy of the soft photon 


k c 


QCD: 


Cancellation of ultraviolet divergence 


Cuv 




Cancellation of infrared divergence 


ClR 




Independence of the cutoff energy of the soft gluon 


k c 



Table 1: List of tests 



3 Calculated processes 

We have calculated the radiative corrections to production processes and decay processes of 
SUSY particles in the framework of the MSSM using GRACE/SUSY-loop. Table OH shows the 
list of processes calculated using GRACE/SUSY-loop. In [3], we have calculated the radiative 
corrections of sfermion-decay processes using the parameter set adopted in [29]. In [4], 



Processes GRACE Preceding studies 



Chargino-pair production (e~ + e + — > \i + Xi) 




12311221 E5] 


Chargino decay (xt ~> two body and xt ^ three body) 


1. 




Neutralino-pair production (e~ + e + — > Xi + X2) 


[2. 


[23JI24] 


Neutralino decay (xl 3 4 — > two body and x° — * three body) 




[26] 


Sfcrmion decay (/ — > two body) 


® 


[2UEH1E9] 


Stop decay (ti->b + xf , t + X? and h ~> b + w+ + X?) 


H[5] 


12SH2S1 


Gluino decay (g —> b + bx, t + t\) 


® 





Table 2: List of processes calculated using GRACE/SUSY-loop 
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however, we use the SPSla' parameter set [3D] for the two-body decay processes of the 
lighter stop (ti) and an original parameter set for the three-body decay process of t\. 



4 Extension of non-linear gauge formalism 

We can extend the NLG functions in the MSSM ([S]) and (jSJ by including bilinear forms of 
sfcrmions with new NLG parameters c's as follows: 



Ft 



w- 



Fz 



{d„ + ieaA„ + igcwPZ li )W +li + i£ w % + 5 H H° + 6 h h° + ikG°)G + 

ij 



(8) 



(<9 M - ieaA^ - igcwPZ^W' 11 - i£vr|(i> + S H H° + S h h° - ikG°)G~ 

ij 

+ E (sr^Si) + ar w + «r*(^)>] - ( 9 ) 

^ + £z 9 -j(v + e H H° + i h h°)G° 



(10) 



where 



= 1,2). 



Processes 



(i,j,M = l,2) 



r, 
n 

e~ 
e~ 
e~ 



"5 

r ; 

-y} 

- T 

- T* 

T* 

3 

Ye 
Ye 
Ye- 



Tk 



T+ + Xl + 



'1 



(1 



Tabic 3: List of tested processes 



We have calculated the cross sections of 
the processes listed in Table |3j in tree-level 
order, and confirmed the NLG invariance of 
the results on the parameters c's in ((HJl, © 
and (|10p. Then we are convinced that the 
extended NLG formalism is valid as a tool 
to test the numerical calculations. 

Figurc[2]shows Fcynman diagrams of the 
process e~ + e + — > + in tree-level 
order (drawn by gracef ig), as an example. 
For this process, coupling constants of the 
vertices, t\ — t-2 — Z and t\ — t-2 — G° (%3) are 
shifted by varying the NLG parameter cJJ > 
but the sum of the amplitudes is invariant. 
Here we set parameters as = 320 GeV, 
Mr 2 = 370 GeV, and total energy as y/s = 
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1000 GeV. The numerical results of the amplitude of each graph and the cross section at 
one point in the phase space are given in Table 2] for (easel) 5iJ=0 and (case2) c[J=1000. 
The total values in two cases agree up to 31 digits, so the NLG invariance of this process is 
confirmed in tree-level order. 

Graph Absolute value of the amplitude 

easel: 1 1.6624897226795816149294531250854308x 10" 3 

2 8.1590691554170053568905511607157850x 10~ 15 

3 1.9404461824036032871287809608655740x 10~ 10 

4 2.0390333192991456823825251175977558xl0~ 9 

5 5.6886720681154757343568160836757535X1Q- 14 

total 1.6624919797058214816061810564292905xlQ- 3 

case2: 1 1.6625481763453484261164965199981440x 10" 3 

2 8.1590691554170053568905511607157850x 10~ 15 

3 1.9404461824036032871287809608655740x 10~ 10 

4 2.0390333192991456823825251175977558X10" 9 

5 5.8453722653531868198152256247220933xlQ- 8 

total 1.6624919797058214816061810564293247X 10~ 3 

Table 4: Tests for the NLG invariance 



5 Summary 

We have developed the program package 
GRACE/SUSY-loop for the EW corrections 
and QCD corrections of the MSSM ampli- 
tudes in one- loop order. Then we have cal- 
culated the radiative corrections to produc- 
tion processes and decay processes of SUSY 
particles in the framework of the MSSM us- 
ing GRACE/SUSY-loop. We have also devel- 
oped a version of GRACE/SUSY-loop for the 
extended NLG formalism by including bi- 
linear forms of sfcrmions with new NLG pa- 
rameters, and have carried out tests for it in 
tree-level order. 
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Figure 2: Feynman diagrams of the process 
e~ + e + — > Ty + T2 in tree-level order. 
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